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ABSTRACT 


The  very  low  frequency  (VLF,  3  to  30  KHz)  part  of  the  radio  frequency 
spectrum  is  characterized  by  low  attenuation  rate,  high  phase  and  frequency 
stability,  and  high  signal  to  noise  ratio.  Consequently,  VLF  radio  propaga¬ 
tion  is  used  for  maty  practical  applications,  e.g.,  frequency  standardization 
clock  synchronization,  and  reliable  long-distance  radio  communications.  Be¬ 
cause  of  the  distinct  advantages  of  VLF  radio  propagation,  the  U.S.  Navy  will 
be  conducting  a  balloon-to-balloon-bome  cross  link  communication  experiment 
to  study  the  characteristics  of  the  VLF  high  altitude  radio  propagation. 

In  this  technical  report,  an  attempt  has  been  made  to  make  theoretical 
computations  of  the  vertical  components  of  the  individual  and  multimode  field 
strengths  as  a  function  of  distance  based  on  mode  theory.  The  variations  of 
various  ionospheric  parameters,  e.g.,  attenuation  rates,  the  height  of  the 
ionospheric  reflection  point,  the  height  gain  factors  for  an  appropriate  com¬ 
bination  of  the  transmitting  and  receiving  antenna  elevations,  along  with  the 
presence  of  the  earth’s  geomagnetic  field,  especially  for  the  East-West 
propagation,  have  been  duly  considered.  The  results  obtained  have  been  pre¬ 
sented  in  tabular  and  graphical  form  and  are  consistent  with  the  values  ob¬ 
tained  by  earlier  workers.  These  field  strength  values  will  be  compared 
against  the  experimental  values  when  the  above-mentioned  ambitious  experiment 
is  carried  out  in  the  Pacific  in  the  late  Summer  of  1983. 
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INTRODUCTION 

Analogous  to  microwave  propagation  along  a  waveguide,  very  low 
frequency  (VLF,  3  to  30  KHz)  radio  propagation  over  long  distance  is 
treated  as  if  the  propagation  takes  place  in  a  terrestrial  spherical 
waveguide  with  boundaries  formed  by  the  conducting  earth's  surface  and 
the  lower  edge  of  the  reflecting  ionosphere.  For  the  frequency  range 
under  consideration,  (15  to  30  KHz,  i.e.,  the  upper  band  of  the  VLF 
spectrum)  the  analysis  of  the  wave  propagation  is  simplified  because 
of  the  limited  number  of  significant  modes  that  need  to  be  considered. 

The  level  of  excitation  of  the  individual  modes  depends  upon  the 
corresponding  excitation  factors  and  the  variation  of  field  strengths 
across  the  waveguide  cross-section  as  given  by  the  so-called  height 
gain  functions.  For  lower  order  modes  (especially  say,  the  first  order 
mode)  the  energy  of  the  wave  is  mostly  concentrated  at  the  base  of  the 
ionosphere.  For  such  a  mode,  the  excitation  factor  is  small,  which 
implies  inefficiency  for  launching  such  a  wave  from  a  ground-based 
point  as  the  excitation  factor  for  such  a  wave  will  not  depend  much 
on  ground  conductivity.  However,  as  it  turns  out,  the  height  gain 
factors  for  such  waves  increase  as  the  height  of  the  transmitting 
antenna  increases.  The  waves  become  ea*.Lk -detached  and  are  called 
earth-detached  modes  (nomenclature  due  to  James  R.  Wait)  or  whispering 
gallery  modes  (nomenclature  due  to  K.  G.  Budden) .  Such  waves,  naturally 
characterized  with  very  low  attenuation  rates,  criss-cross  the  ionosphere 
without  any  intermediate  earth-based  bouncing  point  and  propagate  over 
long  distances  without  losing  much  of  their  energy. 

The  U.S.  Navy  is  planning  a  balloon- to-b alio on  communication  link 
experiment  employing  this  mode  of  propagation.  The  plans  for  this  ex¬ 
periment  have  been  given  elsewhere  (3).  In  this  technical  report,  only 
the  theoretically  computed  values  of  the  vertical  field  components  of 
the  field  strength  radiated  out  by  a  vertically  polarized  antenna  at 
different  distances  for  different  configurations  of  the  transmitting 
and  receiving  antenna  heights  will  be  presented. 
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THEORY 


The  theory  of  very  low  frequency  radio  propagation  has  been  thoroughly 
discussed  in  many  well  known  texts  and  research  reports  (2, 4,5, 6  and  7), 
Several  approximation  theories  can  be  utilized  in  order  to  limit  the 
mathematical  sophistication  within  bounds.  In  this  consideration,  the 
earth-ionosphere  terrestrial  waveguide  is  considered  as  a  thin  spherical 
shell.  Spherical  coordinate  system  (r,0,ij>)  is  used  and  the  boundary 
conditions  at  the  two  interfaces  are  defined  by  their  impedance  conditions, 

1*e*»  E  E 

A  6 

V  H  1  Zi‘  — 

a  a. 


where 


impedance  of  the  conducting  ground 


*  impedance  of  the  conducting  ionosphere 
Eg  «  tangential  electric  field  component 

H<f>  ■  zonal  magnetic  field  component 

Furthermore,  for  special  interest,  discussion  will  be  restricted  to  two 
discrete  frequencies,  19.4  and  27  kHz,  as  these  are  the  two  frequencies 
which  will  be  operational  during  the  balloon-to-balloon-borne  cross  link 
communication  experiment.  Two  levels  of  transmitter  power  will  be  used, one 
kilowatt  and  440  watts.  (For  details,  please  refer  to  (3)  ).  The  ionosphere 
will  be  characterized  with  an  exponential  profile  having  a  conductivity 
•gradient  £»  2  (or  8  *  0.5),  where  is  the  lengt£  over  which  ionospheric 
conductivity,  which  is  a  measure  of  the  ratio  w0/v  or  m/y)t  changes  by  a 

ratio  of  2.71,  where  uQ  is  the  plasma  frequency  and  v  is  the  collisonal 
frequency  and  N  is  the  electron  density  in  the  ionospheric  layer. 

According  to  the  formulations  of  WKB  (or  second  order)  approximations, 
the  integral  form  of  the  modal  equation  or  the  resonance  condition  for 
mode  propagation  of  spherical  waves  can  be  written  as  (2,4,6) 


Rg(Cn)  Ri(c')  exp(-2ik  £  (Cn  +  2z/a)  dz  +  iir/2  ) 


«  exp  (-2imr  )  (1) 

where  R  ( Cn)  and  Rj(Cn)  are  the  Fresnel  reflection  coefficients  of  the 
earth ’s8surf ace  ana  the  ionospheric  surface  for  the  complex  angles  of 


incidences ; 


(1  -  S?  )X/2 


a  -  s'2  )/2 


■:  -v.  *. 


t'. «.■  u tvo v  t^Tvis^Tc-K'K^v's.- *.\»  - «.-•  *r.  *.-  *.-  v  %-  >.-  • 
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The  other  constants  are  defined  as  follows: 

k  *  propagation  constant  =  2it  /X 
a  *  radius  of  the  earth  surface 

zo  *  -  aCn  /2 

As  stated  by  Wait  (4) ,  under  this  approximation,  valid  for  the  frequency 
range  of  our  interest  (20  to  30  kHz),  the  modal  characteristics  are 
independent  of  the  electrical  properties  of  the  earth’s  surface  and  the 
?  modes  are  characterized  as  earth-detached  mode  in  analogy  to  their  counter^ 

parts  in  accoustics  originally  investigated  by  Lord  Rayleigh  at  the  turn 
of  the  present  century. 

*  The  total  phase  given  by  the  integral  form  of  the  modal  equation 

consists  of  a  ”to  and  fro”  path  between  the  ionospheric  reflection  point 
at  z=h  and  a  second  reflecting  point  at  z  =  z  a  height  embraced  by  a 
surface,  called  caustic  surface,  below  which  the  waves  of  lower  order  modes 
becomes  evanescent.  The  phase  shift  at  the  caustic  accounts  for  7r/2, 

The  excitation  factors  for  these  modes  for  the  ground-based  launching  points 
are  very  small. 

Assuming  earth’s  surface  as  a  perfect  elect rtcconductor  and  the 
ionospheric  surface  as  a  perfect  magnetic  conductor  (i.e.,  R  =  1  and 
R.  =  -1),  equation  (1)  reduces  to  8 

h  2  1/2 

-  2ik  /"  (<T  +  2z/a)  dz  +  iir  /2  -  -2ni  ir  +  iir  (2) 

zo 

and  the  solution  is: 

2/3 

cn  “  <  S  (n  "  1M)  )  ~  2h/a 

where  Cq  specifies  the  resonance  condition  for  the  modes. 

The  expression  for  the  amplitude  of  the  ve-tical  component  of  the  electric 
field  strength  at  a  certain  distance  can  be  written  as  (2,6) 


300  P *  x  X  1/2 

Ez  ■  —  <  a  Sin  D/a  >  *  A.  V> 


(3) 


exp  (-  anD) 


volt/meter 


and  the  phase  of  the  nth  mode  component  relative  to  the  antenna  current  is 
given  by  (6) 


E 


tt/4  -  wD/v. 


z,n 


p,n 


(4) 


n 
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where  Ez  *  Modesum  field  strength  amplitude  in  volt/meter 
h  »  ionospheric  height  in  meters 
X  *  c/f  *  wavelength  in  meters 
Pjj  *  radiated  power  in  kilowatts 

D  *  distance  in  meters 

*  nth  mode  excitation  factor 

Gn(z)  -  height  gain  factor  for  the  transmitting  antenna  at 
height  z  meters 

GgCz)  *  height  gain  factor  for  the  receiving  antenna  at 
height  z  meters 

on  *  nth  mode  attenuation  constant  in  dB/megameter 
Ez>n  -  phase  of  the  nth  mode  in  radians 

*n  *  phase  of  the  nth  mode  excitation  factor  in  radians 
b)  *  angular  frequency  of  the  wave  in  radians 
VP  n  “  phase  velocity  of  the  nth  mode  in  meters /second. 

As  noted  earlier,  for  the  frequency  range  under  consideration  (15-30  kBfc), 
very  few  modes  are  important  for  consideration.  Consequently,  for  the 
present  prupose  only  the  first  three  modes  (mode  1  through  3)  will 
be  considered. 
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SINGLE  MODE  FIELD  STRENGTH 


We  can  express  the  field  strength  amplitude  in  dBv/m  from  equation  (3) 
for  the  nth  mode  as  (where  n  is  either  1  or  2  or  3) : 

Ezn  -  20  Log  (300)  -  20  Log  (h)  +  10  Log  PR  +  10  Log  (3X108) 


-10  Log  (f)  -  10  Log  (a  sin  D/a)  -  ctnD 

106 

+  HGF  dBv/m  (5) 

-  20  Log  (106)  +  20  Log  (300)  -  20  Log  (h) 

+10  Log  (Pr)  +  10  Log  (3X108)  -  10  Log  (f) 

-10  Log  (a  sin  D/a)  -  CtnD 

106 

+  HGF  dByv/m  (6) 

For  daytime  condition  (h  *  70  km)  and  for  frequency  f  =*  19.4  kHZ  and 
Pr  »  1  kilowatt,  equation  (6)  reduces  to: 

EZjD  *=  54.3  -  10  Log  (a  sin  D/a)  -  otD 

+  HGF  dByv/m.  kw  (7) 

and  for  nighttime  condition  (h  *  90  km) 

E2jn  "  53.1  -  10  Log  (a  sin  D/a)  -  aD 

+  HGF  dByv/m.  kw  (8) 

Similarly,  for  the  second  frequency,  f  =  27  kHZ  and  for  the  day  and  night¬ 
time  conditions,  the  corresponding  expressions  are: 

Ez>n  =  52.3  -  10  Log  (a  sin  D/a)  -  ou 

+  HGF  dByv/m. kw  (9) 

(f  -  27  kHZ,  h  =  70  km) 


* «  \ 


5 


NADC-83062-40 


Snd 


Ez  n  “  50,9  ”  10  ^(a  sin  D/a)  -  aD  +  HGF  dByv/m.kw  (10) 
(f*27  kHz,  h.=90  km) 


In  equations  (7)  through  (10),  a  and  D  are  expressed  in  megameters  and  a 
is  in  dB/megameter  (Mm) .  HGF  is  the  modified  height  gain  function  for  the 
nth  order  mode  and  is  given  by  (2) 


HGF  «  20  Log  (  ^  Gn(z)  Gn(z)  )  dB 


(ID 


Following  reference' (2)  we  can  write 


where 


4  3/2 


2  t1/2(y  -  t  ) 
n  VJo  n' 


172  <-  3  *H  ) 


(12) 


Gn(z^  *  A 

Ci/4  Ai(tn"z)  «P  <+  |  tn/2  > 

(13) 

and 

Gn(z)  *  A 

i/4  2  3/2 

fcn  Ai(tn~  z)  fcn  } 

(14) 

2/3 

(ka/2)  (2h/a) 


A 

z 


2/3 

(ka/2)  (1^/a) 

2/3 

(ka/2)  (h_,/a) 

2/3  * 


tn  =  (ka/2) 


(-  Cn2) 


C2  »  (  —  (n  -  l/4))2/3 
n 


-  2h/a 


"r 


height  of  the  transmitting  antenna  above  the  ground  surface 
height  of  the  receiving  antenna  above  the  ground  surface 


A^(tn  -  z)  and  A^(t^  -  z)  are  the  Airy  functions  of  arguments 


(tn  -  z)  and  (t  -  z)  as  defined  in  reference  (8) . 


r  — '  _ff  -- 
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HGF  ■  h  (4II4  )1/3  (n-l)~1/3f2^3 
6aC2  4 


X  [Ai(2  (]Ia)2/3f2/3(h-hT)  -  -2.83  (n-l)2/3)J 

4 


a  c 


X  {Ai(2  (IIa)2/3f2/3(h-hR)  -  2.83  (n-l)2/3)] 

4 


a  c 


or 


HGF  *  3.49  X  10-2h  (n-l)“1/3 

4 


X  Ai(3.72  X  10"2(h-hx)  -  2.81  (n-l)2/3) 


X  A  (3.72  X  10~2(h-hR)  -  2.81  (n-l)2/3) 

4 


for  f  *  19.4  KHz 


(17) 


and 


HGF  -  4.35  X  10"2  h(n-l)”1/3  X  Ai(4.64  X  l(f2(h-hT)  -  2.81(n-l)2/3) 

4  4 


X  Ai(4.64  X  10“2  (h-hR)  -  2.81  (n-l)2/3) 

4 


(18) 


for  f  *  27  lBz 


Tabulated  values  of  Airy  functions.  A,  for  different  values  of  their 
arguments  as  obtained  for  various  values  of  hx  and  hR»  the  heights  of  the 
transmitting  and  receiving  antennas»are  listed  in  Table  I.  The  values  have 
been  obtained  from  reference  (8).  The  appropriate  values  of  a(dB#m),  the 
attenuation  constant,  for  different  modes  (mode  1  through  3),  and  for  day 
and  night  time  conditions  are  given  in  Table  II.  These  values  were  obtained 
from  reference  (1). 
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From  the  values  given  in  this  table,  it  is  found  that  the  attenuation 
rate  for  the  first  order  mode  is  less  at  the  lower  frequency  during  the 
daytime  and  increases  during  the  nighttime.  On  the  contrary,  for  the 
second  and  third  order  modes,  the  attenuation  constant  decreases  during 
the  nighttime.  The  trend  is  the  same  for  the  higher  frequency  (27  kHz) . 

Using  the  values  of  Airy  functions  for  various  augments  as  listed  in 
Table  II,  the  corresponding  values  of  modified  height  gain  factors  for 
all  the  modes  (mode  1  through  3)  have  been  calculated  for  the  condition 
when  the  receiving  and  transmitting  antennas  are  at  the  same  elevation. 

The  controlling  parameters  are  the  transmitting/receiving  antenna 
elevation,  the  propagation  condition  for  the  ionosphere,  and  the  trans¬ 
mission  frequency.  The  numerical  values  are  given  in  Table  II  and 
shown  in  graphical  form  in  Figures  1  and  2.  The  calculations  have 
also  been  repeated  for  the  condition  when  the  receiving  antenna  is 
ground-based  and  the  elevation  of  the  transmitting  antenna  is  increased 
above  the  ground.  These  tabulated  values  ate  given  in  Table  III  and 
the  graphical  forms  have  been  represented  in  Figures  3  and  4.  It 
should  be  emphasized  at  this  point  that  these  values  of  the  modified 
height  gain  factors  along  with  the  corresponding  values  of  the  attenuation 
constant  determine  the  values  of  the  vertical  component  of  the  individual 
and  multimode  field  strength  as  a  function  of  distance  and  the  resultant 
model  interference  pattern. 

With  these  calculated  values  of  a  and  HGF's,  the  amplitude  of  the 
vertical  component  of  the  field  strength  for  any  of  the  modes  (mode  1 
through  3),  for  any  frequency  (either  19.4  or  27  kHz)  and  for  any 
propagation  condition  of  the  ionosphere  (h=70  km  for  the  daytime  condition 
and  h=90  km  for  the  nighttime  condition)  can  be  computed  by  using  the 
Equations  7  through  10. 
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Height  (  Ht  =  Hr  ,  f  -  19,4  khz) 


FIGURE  2.  MODIFIED  HEIGHT  GAIN  FUNCTION  VS  ANTENNA  HEIGHT  (H^,  -  HR,  F  -  27  KHZ) 


ANTENNA  HEIGHT  (Hr  OR  Hr)  IN  KM 


Figure  3.  Mddified  Height  Gain  Function  Vs  Antenna  Height.  (Hr  =  0.. ,  f=19.4  khz) 


Figure  4.  Modified  height  Gain  Factor  Vs  /Wtenna  Height  (Hr=0,  f=27  khz) 
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TABLE  I.  VALUE  OF 


AN  AUXILIARY  VARIABLE 


f  (kHz) 


h  Xian) 


a  (dB/Mm) 


^  -3  2/3  i 

(T  =4 ,83x10  h  f  (n--) 


1.56 

6.56 
15.19 


7.2193 

5.443 

4.682 


1.62 

3.20 

7.80 


9.282 

6.998 

6.019 


2.00 

4.54 

10.54 


10.047 

7.575 

6.516 


12.918 

9.740 

8.377 
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TABLE  II.  VALUES  OF  MODIFIED  HEIGHT  GAIN  FACTORS  (frT  =  ^R) 


f 


5 


f  (kHz) 
19.4 


27 


h  (km) 
70 


90 


70 


^T  -(km) 

n 

tn“z 

Ai 

HGF  (dB) 

0 

1 

0.2844 

0.2837 

-  4.72 

2 

-1.4767 

0.4643 

1.39 

3 

-2.9116 

-0.3419 

-  5.24 

10 

1 

-0.0876 

0.3783 

0.28 

2 

-1.8487 

0,3008 

-  6.15 

3 

-3.2836 

-0.4172 

-  1.78 

20 

1 

-0.4596 

0.4674 

3.96 

2 

-2.2207 

0.0962 

-25.96 

3 

-3.6556 

-0.282 

-  8.58 

30 

1 

-0.8316 

0.5266 

6.03 

2 

-2.5927 

-0.1785 

-15.22 

3 

-4.0276 

-0.0703 

-32.71 

0 

1 

1.0284 

0.1353 

-15.4 

2 

-0.7327 

0.5152 

5.38 

3 

-2.1676 

0.096 

-25.12 

10 

1 

0.6564 

0.1973 

-  8.84 

2 

-1.1047 

0.5338 

6.0 

3 

-3.5396 

-0.1123 

-22.39 

20 

1 

0.2844 

0.2834 

-  2.55 

2 

-1.4767 

0.4643 

3.57 

3 

-2.9116 

-0.3419 

-  3.05 

30 

1 

-0.0876 

0.3783 

2.47 

2 

-1.8487 

0.2868 

-  4.80 

3 

-3.2836 

0.4172 

0.40 

0 

1 

0.9284 

0.1468 

-13.29 

2 

-0.8327 

0.5266 

6.45 

3 

-2.2676 

0.0267 

-46.66 

10 

1 

0.4644 

0.2408 

-  4.69 

2 

-1.2967 

0.5123 

5.97 

3 

-2.7316 

-0.2400 

-  8.51 

20 

1 

3.9611X10"4 

0.3550 

2.05 

2 

-1.7607 

0.3408 

-  1.14 

3 

-3.1956 

-0.4175 

1.11 

30 

1 

-0.4636 

0.4674 

6.83 

2 

-2.2247 

0.0961 

-23.10 

3 

-3.6596 

-0.2820 

-  5.71 

0 

1 

1.8564 

0.0431 

-32.40 

2 

0.0953 

0.3318 

0.61 

3 

-1.3396 

0.5123 

6.84 

90 


-  ’.o.*?.;  « 


.Tt  v_v 


^saawaeg»t3iS6g^iigraifc.’s.  #ic*»  n 


* 
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TABLE  II  (CONTINUED) 


f  cte*} 

27 


h  (  kn^ 

90 


^  (kn> 

in 

vz 

Ai  C.^n-. 

10 

1 

1.3924 

0.0834 

2 

-0.3687 

0.4474 

3 

-1.8036 

0.3408 

20 

1 

0.9284 

0.1468 

2 

-0.8327 

0.5266 

3 

-2.2676 

0.0267 

30 

.  1 

0.4644 

0.2408 

2 

-1.2967 

0.5123 

3 

-2.7316 

-0.2400 

HGF  (dB) 

-20.93 

5.80 

-  0.24 

-11.11 

8.63 

-44.48 

2.51 

8.15 

-  6.33 
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TABLE  III. VALDES  OF  MODIFIED  HEIGHT  GAIN  FACTORS  (hR  V  **1) 


(1H. id  h  ^m.)  Tia  Ckm> 


HGF  (HB) 


1 

-  4.72 

2 

1.39 

3 

-  5.24 

1 

-  2.22 

2 

-  2.38 

3 

-  3.51 

1 

-  0.38 

2 

-12.27 

3 

-  6.91 

1 

0.66 

2 

-23.42 

3 

-18.98 

1 

-15.40 

2 

5.38 

3 

-25.12 

1 

-12.12 

2 

5.69 

3 

1 

-  8.97 

2 

4.48 

3 

1 

-  6.47 

2 

0.29 

3 

-12.35 

1 

-13.29 

2 

6.45 

3 

-46.66 

1 

-  8.99 

2 

6.21 

3 

1 

-  5.62 

2 

2.67 

3 

1 

-  3.23 

2 

-  8.33 

3 
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TABLE  I^I(  CONTINUED) 


(Khz)  h  (tart) 

•  27  90 


10 

20 

30 


_n 

HGF  (dB) 

1 

-32.40 

2 

0.61 

3 

6.84 

1 

-26.66 

2 

3.20 

3 

3.30 

1  • 

-21.75 

2 

4.62 

3 

-18.82 

1 

-17.45 

2 

4.38 

3 
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MULTIMODE  FIELD  STRENGTHS 


At  long  distances  the  field  strengths  due  to  second  and  third  order 
mode  are  equally  important.  Hence,  total  field  strengths  at  large  distances 
are  the  vector  sum  of  individual  mode  field  components.  To  calculate  such 
resultant  field  strength,  a  knowledge  of  the  phase  of  the  individual  mode  at 
a  certain  distance  is  essential. 

The  phase  of  the  vertical  component  of  the  nth  mode  (n  =  l,or  2  or  3) 
field  strength  with  respect  to  the  transmitting  antenna  current  is  given  by 
equation  (4)  as 


TT 


=  —  -  COD/  V 


■■z,n 


p,n 


+  <f>. 


n 


TT 

4 


Dc 


X 


v 

p,n 


+ 


where  c  is  the  speed  of  light  in  free  space  and  the  other  terms  have  their 

previously  defined  meanings.  The  values  of  c/v  and  (j> .  ,  as  obtained 

An 

from  reference  (l),are  given  in  Table  III.  With  the  values  of  c/v  and 

p  ,n 

<f>A  t  values  of  for  any  mode,  the  frequency,  distance  and  ionospheric 

conditions  can  be  calculated  from  the  equation  cited  above.  The  values  of 

^Ez,n  thus  calculated  for  different  modes  and  frequencies  are  listed  in 

Table  IV.  From  this  table,  it  is  also  found  that  the  difference  in  phase 
between  first  and  second  order  modes  for  f  =  19.4  kHz  is 


A 


5.01  radian/Mm  during  daytime  (h  =70  km) 
3.20  radian/Mm  during  nighttime  (h  =90  km) 


and  for  f  =  27  kHz 

A  <|>i  2  =  4.22  radian/Mm  during  daytime  (  h=70  km) 

2.33  radian/Mm  during  nighttime  (h=90  km) 

Furthermore,  with  known  values  of  „  ,  the  amplitude  of  the  vertical 

tz,n 

component  of  the  multimode  field  strength  is  calculated  as  follows: 


E 


(  (Ezl  cos  $Ezi  +  Ez2  cos  (j,^  +  Ez3  cos  <J>ez3  )2 

<  Ezl  sin  <J>Ezl  +Ez2sin  hz2  +Ez3  Sin  ^Ez3  )2  } 
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TABLE 

IV. 

*  VALUES  FOR  <J>‘ 

n>  vpn 
c 

AND  <j>n 

f  (kH2) 

h  (km) 

a 

In 

Vpn 

jfcn 

C 

19.4 

70 

l 

7.50 

0.9984 

0. 916-0. 407D 

2 

1.40 

1.011 

0. 809-0. 402D 

3 

0.78 

1.042 

0. 799-0. 390D 

90 

1 

13.12 

0.996 

1. 014-0. 408D 

2 

1.82 

1.0038 

0. 817-0. 405D 

3 

0.66 

1.02 

0. 797-0. 398D 

27 

70 

1 

14.2 

0.997 

1. 033-0. 567D 

2 

1.6 

1.0035 

0. 813-0. 563D 

3 

0.42 

1.017 

0. 792-0. 556D 

90 

1 

25.5 

0.9944 

1.23-0.568D 

2 

6.5 

0.9994 

0. 898-0. 566D 

3 

-0.03 

1.0065 

0. 785-0. 562D 

is  calculated  from 

$  *  "I  “  *f  $  n  rad 
4  V 


=  7T  -  2_  D  c_  •+*  4>  n 
4  A  Vpn 


where  A,  D  inkn»  ,  A  •  300,  f  in  K^z. 

f 
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and  the  corresponding  phase  is  calculated  from 


=  tan 


Ezl  sin  0Ezl  +  Ez2  sin  0Ez2  +  Ez3  sin  0Ez3  ^ 
EZ1  cos  ^Ezl  +  Ez2  C0S  0Ez2  +  Ez3  cos  0Ez3 


Ez3  sin  0Ez3 
■F3  cos  0Ez3 
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DISCUSSION 


A  microcomputer  (Tektronix  4054  with  accompanying  Hard  Copy  Unit  4631) 
was  used  for  all  the  computational  works  reported  in  this  report.  A  computer 
program  written  in  BASIC  and  provided  in  Appendix  A  has  been  used  to  compute 
the  values  of  E2n,  the  vertical  component  of  the  individual  mode  electric 
field  strength  as  a  function  of  distance  based  on  any  of  the  Equations  7 
through  10.  The  sample  program  shown  in  Appendix  B  was  used  for  the  computation 
of  the  values  of  Ez,  the  multimode  field  strength  based  on  Equation  19  and 
provides  the  output  in  graphical  and  tabular  forms.  These  graphical  plots 
have  been  represented  in  Figures  5  through  20.  For  all  these  plots,  the 
transmitting  and  receiving  antennas  are  assumed  to  be  at  the  same  elevation 
and  the  controlling  parameters  use  either  the  frequency,  the  height  of  the 
ionospheric  reflection  point,  or  the  elevation  of  the  transmitting/receiving 
antennas.  Distances  are  varied  up  to  3000  km  and  the  elevation  of  the  trans¬ 
mitting/receiving  antennas  up  to  30  km  in  order  to  conform  to  the  plans  of 
the  experiment.  Figures  21  through  32  illustrate  similar  results  when  the 
receiving  antenna  is  ground-based  and  the  transmitting  antenna  elevation 
is  varied  up  to  30  km  in  steps  of  10  km. 

The  signature  of  strong  interference  phenomenon  due  to  individual  modes 
on  the  modesum  field  strength  is  apparent  and  it  is  indicative  of  the  fact 
that  at  many  locations,  although  the  individual  mode  field  strength  is  strong, 
the  modesum  field  strength  is  relatively  weak  due  to  modal  interference. 

At  low  frequency  (19.4  kHz)  and  during  daytime,  the  interference  pattern 
begins  to  appear  at  short  distances.  For  the  ground-based  transmitting  and 
receiving  antennas,  the  modesum  field  strength  curve  displays  strong  inter¬ 
ference  pattern  with  bit  nulls  occurring  at  600  km  and  2000  km  as  can  be 
seen  in  Figure  5.  The  interference  phenomenon  becomes  weaker  as  the  antenna 
elevation  is  increased  and  at  an  antenna  height  around  30  km  it  disappears 
completely  as  is  seen  in  Figure  8.  This  variation  of  interference  phenomenon 
can  be  attributed  to  the  fact  that  at  higher  altitude,  the  height  gain  factor 
for  the  first  order  mode  increases  and  for  the  second  and  third  order  mode 
it  decreases.  Furthermore,  the  attempt ion  constant  for  the  first  order 
mode  is  very  low  (1.5  dB/Mm)  whereas  for  the  second  and  third  order  modes 
attenuation  constant  is  fairly  high  (5.6  dB/Mm  and  15.2  dB/Mm).  All  these 
factors  combine  together  to  make  the  field  strength  curve  for  the  first  order 
mode  predominant  at  large  distances  thereby  making  the  interference  pattern 
very  weak. 

Interference  patterns  during  nighttime  are  illustrated  in  Figures  9 
through  12.  It  is  seen  in  these  figures  that  at  nighttime,  the  interference 
phenomenon  is  very  weak  when  the  antenna  elevation  is  low  and  gets  stronger 
as  the  elevation  is  increased  and  is  strongest  when  the  antenna  height  is 
30  km.  A  prominent  null  occurs  at  1000  km  as  is  seen  in  Figure  12.  Referring 
to  Figure  1,  it  is  found  that  at  nighttime  the  attenuation  rate  for  the 
first  order  mo  .e  is  higher  than  during  the  daytime  whereas  for  the  second 
and  third  order  modes,  this  rate  is  lower.  The  height  gain  factor  for  the 
second  order  mode  dees  not  fall  off  fast  with  antenna  elevation  and  for  the 
third  order  mode  it  rather  increases.  The  net  effect  of  these  factors  is 
to  make  the  second  and  third  order  mode  field  strengths  very  strong  and 
results  in  a  strong  interference  pattern  for  the  multimode  field  strength 
curve. 
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At  the  higher  frequency  (27  kHz)  and  during  the  daytime  propagation 
condition,  a  reverse  condition  for  the  interference  is  displayed  by  the  multi- 
mode  field  strength  curve  as  is  shown  in  Figures  13  through  16.  At  lower 
antenna  elevation,  the  second  order  mode  curve  is  weak  and  the  third  order 
mode  curve  is  weakest  making  the  multimode  interference  phenomenon  less 
prominent.  As  the  antenna  elevation  increases,  second  and  third  order  mode 
curves  become  stronger  and  a  strong  interference  pattern  is  observed  with 
a  deep  null  occurring  at  850  km  when  the  transmitting  and  receiving  antenna 
elevation  is  20  km  as  is  illustrated  in  Figure  15.  As  Figure  16  shows,  above 
this  height  the  interference  phenomenon  becomes  weak. 

Figures  17  through  20  show  the  interference  pattern  for  the  higher  fre¬ 
quency  (27  kHz)  during  the  nighttime  condition.  For  the  ground-based  antennas 
strong  modal  interference  accompanied  by  two  deep  nulls  occurs.  The  intensity 
of  interference  decreases  as  the  transmitting/receiving  antenna  elevation  is 
increased . 

The  interference  phenomenon  is  also  present  in  the  multimode  field  strength 
curve  as  a  function  of  distance  both  during  the  daytime  and  nighttime  when 
the  receiving  antenna  is  ground-based  and  the  transmitting  antenna  is  elevated. 
The  interference  pattern  changes  from  weak  to  strong  and  vice  versa  when  this 
elevation  is  increased  up  to  30  km.  As  an  example  for  the  lower  frequency, 
it  is  seen  in  Figure  21  that  at  closer  distances  and  at  lower  antenna  elevations 
the  height  gain  factors  for  the  second  and  third  order  mode  are  fairly  high 
making  their  corresponding  field  strengths  strong.  Interference  phenomenon 
begins  to  appear  at  short  distances,  and  a  strong  null  occurs  at  600  km. 

At  higher  transmitting  antenna  elevation,  the  interference  pattern  fades 
away  as  the  second  and  third  order  mode  fields  become  less  predominant. 

At  nighttime,  the  interference  phenomenon  is  not  so  strong. 

At  higher  frequency  and  during  the  daytime  condition,  the  interference 
phenomenon  is  strong  for  the  intermediate  elevation  of  the  transmitting 
antenna,  the  receiving  antenna  being  ground-based  and  is  characterized  by  the 
occurence  of  two  strong  nulls,  one  at  800  km  and  the  other  at  2400  km  as  are 
shown  in  Figure  24.  During  nighttime,  a  strong  interference  phenomenon  occurs 
only  when  the  transmitting  antenna  elevation  is  around  10  km.  At  higher 
elevations,  the  second  and  third  order  mode  become  less  significant  and 
the  interference  pattern  disappears  -  making  the  multimode  field  strength 
curve  equivalent  to  the  first  order  mode  field  strength  curve. 


FIGURE  5  .  FIELD  STRENGTH  VS  DISTANCE 

H=70  KM, HT=HR=0  KM,  F=  19.4  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
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FIGURE  8  .  FIELD  STRENGTH  VS  DISTANCE 
H=70  KM,  HT  =HR=30  KM, Fa  19.4  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/ttETER . tAW 
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FIGURE  9  .  FIELD  STRENGTH  VS  DISTANCE 
H=90  KM, HT=HR=00  Kf1,F=  19.4  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KV 


FIGURE  15 .  FIELD  STRENGTH  VS  DISTANCE 
H=70  KM,HT=HR=20  KM,  F=  27.0  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KV 


FIGURE  16  .  FIELD  STRENGTH  VS  DISTANCE 
H=70  KM, HT*HR=30  KM, F=  27.0  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KW 


FIGURE  17  .  FIELO  STRENGTH  VS  DISTANCE 
H=90  Kf1,HTsHR=00  Kf 1,F=  27.0  KHZ 
X-AXIS  DISTANCE  IN  MECArtETERS 
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FIGURE  18  .  FIELD  STRENGTH  VS  DISTANCE 
H=90  KM, HT=HR= 1 0  KM,F=  27.0  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KW 


NADC-83062-40 


o  X1 
CM  ,| 

zz 

UJsCi 

tr  ' 
3S: 
O  CD 
»->  II 
U.X 


FIGURE  21  •  FIELD  STRENGTH  VS  DISTANCE 
H=70  M1,HR*0  Ml,  HT=10  M1,F=19.4  KHZ 
X-AXIS  DISTANCE  IN  flEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KW 
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IGURE  23.  FIELO  STRENGTH  VS  DISTANCE 
ss70  KM,  HRss0  KM,  HT=30  KM, F* 19. 4  KHZ 
-AXIS  DISTANCE  IN  MEGAMETERS 
-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KW 


FIGURE  25  .  FIELD  STRENGTH  VS  DISTANCE 
H=90  KM,HRs0  KM,  HT=20  KM,F=19.4  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KV 


FIGURE  27  .  FIELD  STRENGTH  VS  DISTANCE 
H*70  KM,HR=0  KM,  HT=:!0  KM,F*27.0  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER.KV 
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FIGURE  28  .  FIELD  STRENGTH  VS  DISTANCE 
H=70  KM,HR=0  KM,  HT=20  KM,F=*27.0  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER.KV 
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FIGURE  30 .  FIELD  STRENGTH  VS  DISTANCE 
H=90  KM, HR=0  KM,  HT=I0  KM,F*27.0  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KW 


FIGURE  31.  FIELD  STRENGTH  VS  DISTANCE 
H= 90  KM,HR=0  KM.  HT=20  KM,F*27.0  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  OB/uV/METER . KW 
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INFLUENCE  OF  GEOMAGNETIC  FIELD  ON  MODESUM 
AND  INDIVIDUAL  MODE  FIELD  STRENGTHS 


In  all  the  previous  calculations,  it  has  been  tacitly  assumed  that 
the  ionosphere  is  isotropic  and  is  characterized  by  an  exponential  profile 
for  electronic  density  and  collisional  frequency  variations.  This  assumption 
neglects  the  presence  of  the  earth’s  geomagnetic  field.  When  the  presence 
of  this  field  is  taken  into  consideration,  the  ionosphere  becomes  anisotropic 
and  the  characteristics  of  VLF  radio  propagation  change,  though  in  most  cases 
these  changes  are  very  small.  The  presence  of  the  geomagnetic  field  produces 
changes  in  the  values  of  Individual  mode  and  modesum  field  strengths  through 
the  changes  in  the  values  of  phase  velocity,  phase  of  the  excitation  factor 
and  attenuation  constant,  especially  when  the  propagation  takes  place  along 
the  East-to-West  or  West-to-East  directions.  It  has  been  experimentally 
observed  that  the  field  strengths  for  the  East-to-West  propagation  is  weaker 
than  for  the  West-to-East  propagation,  in  clear  violation  of  the  principle 
of  reciprocity  for  waveguide  mode  of  propagation  (9) . 

The  plans  for  the  balloon-to-balloon-borne  communication  experiment 
calls  for  one  of  the  links  to  have  East-to-West  propagation  (Kwajalein,  the 
transmitting  site  to  Guam,  the  receiving  site) .  In  order  to  make  an  estimate 
of  the  influence  of  the  geomagnetic  field  on  this  communication  link,  some 
of  the  calculations  -  especially  for  the  predominant  mode  (the  first  order 
mode)  -  have  been  repeated  for  an  anisotropic  ionosphere.  The  calculations 
are  applicable  for  a  field  station  with  a  dip  angle  of  7°,  which  is  the  dip 
angle  midway  between  Kwajalein-Guam  link  (data  courtesy  of  the  Geomagnetic 
Division,  World  Data  Center  A,  Boulder,  Colorado).  The  values  of  the  various 
parameters  for  an  anisotropic  ionosphere  are  listed  in  Table  V.  Field  strengths 
for  the  first  order  mode  at  19.5  kHz  as  a  function  of  distance  for  both  the 
isotropic  and  anisotropic  ionosphere  have  been  presented  in  graphical  forms 
in  Figures  33  through  36.  Results  for  the  daytime  condition  with  the  trans¬ 
mitting  and  receiving  antennas  having  the  same  elevations  (either  ground- 
based  or  30  km)  are  shown  in  Figures  33  and  34.  Figures  35  and  36  represent 
similar  results  for  the  nighttime  condition.  In  all  these  figures,  the 
top  curve  is  for  the  isotropic  ionosphere  and  the  bottom  curve  for  the  aniso¬ 
tropic  ionosphere.  The  corresponding  numerical  values  of  the  field  strengths 
and  the  differences  of  the  field  strengths  are  listed  in  tabular  forms  in 
Table  VI  through  XVII.  In  all  these  analyses,  it  is  found  that  the  field 
strength  differences  occur  with  a  maximum  value  of  about  3  dB  at  -a  distance 
3  Mm.  Thus,  the  anisotropy  of  the  ionosphere,  though  important,  is  not 
very  significant. 
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TABLE  VI.  VALUES  OF  FIRST  ORDER  MODE  FIELD  STRENGTH  VS  DISTANCE  FOR  ISOTROPIC  IONOSPHERE 
(  h™  =  h  =  0,  h  =  70  km,  f  =*  19.4  kHz  ) 
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TABLE  VII.  VALUES  OF  FIRST  ORDER  MODE  FIELD  STRENGTH  VS  DISTANCE  FOR  ANISOTROPIC  IONOSPHERE 
(  hT  =  h»  0,  h  =  70  km  ,  f  ■  19.4  kHz  ) 


0.Z2  2  0.04  3  0.06  4  0.09 
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FIELD  STRENGTH  DIFFERENCE  (DB . uV/METER . KV)  VS  DISTANCE  (NUMBER  X  25  KM) 
F=  1 9 . 4  KHZ,  H=70  KM,  HT=HR=0  KM,  n=1 

TABLE  VIII. 


FIGURE  34  ..FIELD  STRENGTH  VS  DISTANCE 
H=70  KM,HR=30KM,  HT=30  KM,F=19.4  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER. KV 


TABLE  IX.  VALUES  OF  FIRST  ORDER  MODE  FIELD  STRENGTH  VS  DISANCE  FOR  ISOTROPIC  IONOSPHERE 
(  hT  m  h„  ■  30  km,  h  =  70  km,  f“  19,4  kHz  ) 
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TABLE  X.  VALUES  OF  FIRST  ORDER  MODE  FIELD  STRENGTH  VS  DISTANCE  FOR  ANISOTROPIC  IONOSPHERE 
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TABLE  XI.  VALUES  OF  THE  DIFFERENCE  BETWEEN  FIRST  ORDER  MODE  FIELD  STRENGTH  VS  DISTANCE 
FOR  ISOTROPIC  AND  ANISOTROPIC  IONOSPHERE  . 

(  hm  =  h„  «  30  km,  h  «  70  km,  f  <=>  19.4  kHz  ) 


FIGURE  35  .  FIELD  STRENGTH  VS  DISTANCE 
H=90  KM, HRs  0KM,  HT=  0  KM,F=19.4  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  OB/uV/ftETER.KW 
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TABLE  XIX.  VALUES  OF  FIRST  ORDER  MODE  FIELD  STRENGTH  VS  DISTANCE  FOR  ISOTROPIC  IONOSPHERE 
(  h_  »  hD  -  0,  h  =  90  km,  f  *  19.4  kHz  ) 
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TABLE  XIII,  VALUES  OF  FIRST  ORDER  MODE  FIELD  STRENGTH  VS  DISTANCE  FOR  ANISOTROPIC  IONOSPHERE 
(  h  =  hD  =  0,  h=  90  km,  f  -  19.4  kHz  ) 
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FIGURE  36  .  FIELD  STRENGTH  VS  DISTANCE 
H=90  KM, HR=30KM,  HT=30  KM,F«19.4  KHZ 
X-AXIS  DISTANCE  IN  MEGAMETERS 
Y-AXIS  FIELD  STRENGTH  IN  DB/uV/METER . KW 
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TABLE  XV.  VALUES  OF  THE  FIRST  ORDER  MODE  FIELD  STRENGTH  VS  DISTANCE  FOR  ISOTROPIC  IONOSPHERE 
(  hT  *  h  =  30  km,  h  «  90  km,  f  »  19.4  kHz  ) 
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FOR  ISOTROPIC  AND  ANISOTROPIC  IONOSPHERE. 

(  hm  =  h_  =  30  km,  h  =  90  km,  f  =■  19.4  kHz 
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NADC -83062-40 

CONCLUSION 


A  comprehensive  computational  analysis  of  the  vertical  component  of 
the  electric  field  strength  for  individual  and  multimode  propagation  of 
electromagnetic  waves  at  very  low  frequencies  (19.4  KHz  and  27  KHz)  as  a 
function  of  distance  for  different  ionospheric  conditions  and  transmitting 
and  receiving  antenna  elevation  configurations  has  been  completed.  The 
analysis  proves  that  the  field  strength  value  for  the  individual  modes 
is  dependent  on  the  attenuation  constant  and  the  modified  height  gain 
factor  which  have  been  proved  to  be  dependent  on  the  transmitting  and 
receiving  antenna  elevation  configurations.  Consequently,  at  moderately 
large  distances,  fields  of  higher  order  modes  beomce  as  significant  as 
that  of  the  first  order  mode. 

The  multimode  field  -strength  curve  displays  strong  interference 
phenomenon.  The  intensity  of  this  interference  at  a  certain  distance 
is  dependent  upon  the  transmitting  frequency,  ionospheric  condition  and 
the  transmitting  and  receiving  antenna  height  configurations.  When  the 
model  interference  is  strong,  the  interference  pattern  displays  the 
occurrence  of  deep  nulls  which  sometimes  occur  at  distances  800  Km  and 
2400  Km  especially  at  nighttime  when  both  the  antennas  are  ground-based 
and  the  frequency  of  transmission  is  27  KHz. 


NADC-83062-40 


REFERENCES 


C.  E.  Brookes,  Jr.,  J.  H.  McCable,  and  F.  J.  Rhoads,  Theoretical  VLF 
Multimode  Propagation  Predictions  NRL  Report  6663,  December,  1967 

Peter  R.  Bannister,  VLF/LF/MF  Whispering  Gallery  Propagation  Studies, 
NUSC  Technical  Report  6503,  September,  1981 

E.  Wykes,  Test  Plan  -  Strategic  Balloon-Borne  Communication  Experiment, 
NAVAIRDEVCEN  Work  Unit  Number  CY  930,  February,  1983 

James  R.  Walt,  Electromagnetic  Waves  in  Stratified  Media,  Pergamon 
Press,  Oxford,  1970 

James  R.  Wait,  Concise  Theory  of  Radio  Transmission  in  the  Earth- 
Ionosphere  Waveguide,  Reviews  of  Geophys  and  Space  Phys,  16  No.  3, 
320-326,  1978  ‘  ” 

A.  D.  Watt,  VLF  Radio  Engineering,  Pergamon  Press,  Oxford,  196/ 

Janis  Galegs,  Terrestrical  Propagation  of  Long  Electromagnetic 
Waves,  Pergamon  Press,  Oxford,  1972 

M.  Abramowitz  and  I.  A.  Stegum,  Handbook  of  Mathematical  Functions, 
National  Bureau  of  Standards,  Washington,  DC,  1964 

Ya  L.  Al'Pest  and  D.  S.  Fligel,  Propagation  of  ELF  and  VLF  Waves  Near 
The  Earth,  Consultants  Bureau,  New  York,  1970 


APPENDIX  B.  PROGRAM  FOR  MULTIMODE  FIELD  STRENGTH  COMPUTATION 
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